Quantum Versus Mean Field Behavior of Normal Modes of a Bose-Einstein 

Condensate in a Magnetic Trap 
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Quantum evolution of a collective mode of a Bose-Einstein condensate containing a finite number 
N of particles shows the phenomena of collapses and revivals. The characteristic collapse time de- 
pends on the scattering length, the initial amplitude of the mode and N. The corresponding time 
values have been derived analytically under certain approximation and numerically for the parabolic 
atomic trap. The revival of the mode at time of several seconds, as a direct evidence of the effect, 
can occur, if the normal component is significantly suppressed. We also discuss alternative means 
to verify the proposed mechanism. 
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The progress in trapping the low density alkaline 
gases and cooling them below the temperature Tbec of 
the Bose-Einstein condensation has initiated a search for 
features of the atomic coherent collective behavior spe- 
cific for the trapped condensate. The theoretical predic- 
tions made for the normal frequencies of the conden- 
sate in the traps have been successfully confirmed in the 
experiments ]8|,^| employing the resonant modulations of 
the trapping potential. 

It is worth noting that the calculations of the normal 
frequencies are based on the Gross-Pitaevskii (GP) 
equation strictly valid in the limit of a number of particles 
N — > oo, with the average density p being fixed flO|| . The 
GP equation describes the condensate in the mean field 
(MF) approximation. In this sense the condensate wave 
function $(x, t) obeying the time dependent GP equation 
is viewed as a classical field [[l0| . In the atomic traps 0-^] 
neither N can be obviously considered as infinite nor p 
is fixed (p ~ N for N < 200 — 500 and p ~ N 2 / b for 
N > 10 4 |pT| ). Therefore, it is important to understand 
how the actual quantum evolution of the collective mode 
deviates, if any, from that predicted by the mean field 
approach for finite N. 

In this paper we will show that, while the normal fre- 
quencies and short time evolution are determined cor- 
rectly by the GP equation whenever N is only a few times 
larger than 1 , the quantum behavior of the normal modes 
can deviate from the mean field one at times much longer 
than the time scale set by the normal frequencies. The 
quantum behavior is characterized by dephasing, that is 
by the phenomena of the collapses and the revivals of 
the normal mode amplitude at zero temperature. Below 
we will derive the corresponding quantum solution in the 
approximation neglecting the exchange effects. The nu- 
merical solution taking into account the full many body 
Hamiltonian for N = 300 in the approximation allowing 
N atoms to occupy two single particle levels only is pre- 
sented also. It turns out that both approaches yield close 
values for the collapse and the revival times, respectively, 
with the collapse time being of the same order of magni- 
tude as the relaxation time reported experimentally . 



The quantum collapses and revivals were first analyzed 
in Ref. |0| for a single atom in a resonant cavity. Very 
recently such phenomena have been observed experimen- 
tally 1 3-3 1 . As was suggested in Refs. JlJJl^], the phase 
of the condensate containing finite numbers of particles 
in their ground state should also undergo the same phe- 
nomena. In this case, however, the one-particle density 
matrix remains insensitive to such fluctuations. Conse- 
quently, special arrangements of at least two condensates 
are required in order to detect this effect jl4|,|l6|]. 

In contrast, for the condensate collective mode the sit- 
uation is less restrictive. The collapses and revivals of 
the collective mode of the condensate with fixed N can 
be detected by the one-particle density matrix at zero 
temperature. A physical reason for this effect is the fol- 
lowing: For fixed finite N a collective mode is represented 
by a multiple mixture of the eigenstates rather than by 
a single eigenstate of the TV-body Hamiltonian. In gen- 
eral, the inter particle interaction makes the eigenener- 
gies dependent nonlinearly on the quantum numbers of 
the Fock space. Under such circumstances the quantum 
beats develop resulting in collapses and revivals [jl7| . For 
N — > co , this mixture approaches the coherent state and 
can become a true eigenstate of the Hamiltonian. Ac- 
cordingly, in such a limit the beats should vanish. It is 
worth noting that in the atomic traps the existence of 
this limit is nontrivial because of the dependence of p on 
N. 

Within the Hartree approach excitations of a Bose- 
Einstein condensate can be described by the GP equation 
i$ = for the classical field $(x,t) 

which represents the Hartree many body wave function 
|$) for N bosons. Here = is the GP energy 

functional obtained from the second quantized Hamilto- 



(1) 



by replacing the Bose operator obeying the usual Bosc 
commutation rule, by the classical field $ jlC| ]. As usual, 
in Eq.(l) Hi is the single particle Hamiltonian; uq de- 
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notes the interaction constant (we express uq = ina/M 
by means of the positive scattering length a and the 
atomic mass M in units where K = 1); p is the chem- 
ical potential. In fact, the replacement of the operator 
^ by $ is equivalent to projecting of the full many body 
wave function |\&) to the Hartree (MF) ansatz |U| 

m = E^|m i £d(tx7 fc (t) - 1, (2) 

with Ck(t) being c-numbers and |0) standing for the vac- 
uum state. In (2) a\, are creation and annihilation op- 
erators, respectively, for a particle at the single particle 
state ipf. (x) . For the most general representations of the 
classical field $(x, t) — J2k CfcWV'fcM and the Heisen- 
berg operator \&(x, t) — J2k a k(t)^k(^-), the replacement 
procedure and calculating the mean ($|ff |$) for the state 
(2) yield the same GP functional E^> from the Hamil- 
tonian (1) when TV is only a few times larger than 1. 

In what follows we will compare predictions for the 
time evolution of the normal modes provided by the mean 
field approach with the full quantum mechanical calcu- 
lations which do not rely on the projection to the states 
(2) . A good qualitative description of the quantum evo- 
lution is provided by the following two-level model: N 
bosons can occupy two one particle states, one given by 
the function V'o(x) minimizing the Hartree energy in the 
ground state and the other given by the function ipi(x) 
which is orthogonal to ipo(x). Both functions are taken 
orthonormal J dxip*ipj = Sij. Keeping i,j = 0,1 only 
restricts the space of states available for each particle. In 
this manner, we neglect the processes of interaction of the 
collective mode under consideration, which corresponds 
to transitions between states ?/>o( x ) and ^i(x), with the 
rest of other modes. Note that the interaction with other 
modes, ascribed to transitions to states ijjk with k > 1, 
is sensitive to the square of the mode amplitude, while 
the rate of the dephasing, as it will be shown, depends 
linearly on the amplitude. Therefore, for the small am- 
plitude the dephasing effect will dominate. 

The choice of ip® (x) is dictated by the requirement that 
the Hartree variational energy of the ground state is mini- 
mal for Co = 1 (and Ck = for k > 0) in (2). Taking into 
account Eqs.(l),(2) and minimizing E^> with respect to 
$ = %/jo, one obtains the GP equation 

(H 1 -v)iP a + u N\^o\ 2 ipo = 0- (3) 

The choice of -01 ( x ), however, depends on the char- 
acteristic of the collective mode being studied as well 
as the magnitude of its frequency p9| . In what fol- 
lows we will employ the Schrodinger representation re- 
lying on the time- independent operator Substituting 
\& = aoV'o(x) + ai'0i(x) in Eq.(l), one obtains the second 
quantized Hamiltonian (1) rewritten as 

H = {s Q - p)%ao + (en - n)a\ai + [(e w + 

(4) 

v(t))a\a + h.c] + \ J2i,j,k,l u ijH a i a ] a kau 



where u ljk i = uo J dx.ip*ip*^ k ipi (i,j,k,l = 0,1); e i3 
j dxipi H±ipj . We have also introduced the term ~ v(t) 
which models the effect of the modulation of the trapping 
potential employed in Refs. HQ to excite the condensate 
normal modes. In what follows we assume that ipi are 
chosen in a manner which yields Uijki real. This, in turn, 
implies the condition u^ki = UkUj ■ 

If the total number of particles N is fixed, the Hamil- 
tonian (4) becomes a JV + 1 by A + 1 Hermitian ma- 
trix H n m = (N — n, n\H\m, N — to) in the Fock space 
\n,N — n) where n — 0,1. ..AT denotes the number of 
atoms at the level 1, with N — n atoms remaining on the 
level 0. Correspondingly, the Schrodinger equation for 
the many body wave function in the Fock representation 
I*) =En= B n (t)\n,N-n), EtoW = 1 is 

N 

iB n — S H nm B m , (5) 

m=0 

where explicitly the nonzero elements (H* m = 
H m ,n) are = en + gin 2 / 2 ; ff„,„_i = 

y/n(N - n + l)[v(t) + (nino - u woo )(n - 1)], a nd 
H n , n -2 = unoovM" -l){N-n + l)(N - n + 2)/2. 
Here the notations g\ = uoooo + umi — 4uioio, e = 
En - £oo + (uoooo - Uim)/2 + (2uioio - "oooo) N are intro- 
duced, and the overall constant in H n ^ n is omitted to set 
i?o,o = 0. Note that the term gin 2 /2 in H n ^ n is due to 
the interparticle interaction. This describes the effect of 
renormalizing of the ideal gas spectrum en in the lowest 
order with respect to the scattering length a. We have 
also made use of the relation £io+ Aitiooo — obtained by 
means of multiplying Eq.(3) by ipi and integrating over 
x, given the condition J dx^ipo = 0. The condensate in 
this model is associated with the ground eigenstate of the 
Hamiltonian matrix in Eq.(5) for v = 0. The one particle 
density matrix p(x,x',f) = 'Eij Pij(t)ip*(x)ipj(x.'), with 

Pij(t) = (*|a|aj-|*) (Pij(t) = p*ji{t)), is given explicitly 
by 

p 00 (t) = N — p n (t), Pll (t) = En=o n\B n {t)\\ 



Pio(*) = E„=o y/(n + l)(N-n)B n+1 (t)*B n (t). 

Note that the corresponding MF expression for the 
density matrix is p\f\t) = ($|at aj -|$) = NC*(t)Cj(t), 
where the Hartree ansatz (2) has been employed and 
Ci{t) must be found from the GP equation iCi = 
SE^/SC* in the restricted basis tpo, ip\ 

UijkiCjCkCi. (7) 

j jkl 

Here i, j, k, I = 0, 1; also we have introduced the notations 
v w = Vqi = v(t). This system describes the evolution of 
the Hartree state (2) |Q which is expected to be very 
close to the actual state for the case N — > oo and p fixed 
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Jj"o| ]. For finite and large N one can expect that the 
classical and the quantum solutions do not differ much 
at short times. Therefore, it is natural to assume that 
the driving term v(t) in Eq.(7) prepares the state (2) (at 
t = 0) with some Cfc(0) and then is removed so that 
the quantum evolution at t = starts from the state 
(2), with v(t) — for t > 0. This provides the initial 
condition |\P(f = 0)) = |$(t = 0)). Comparing (2) with 
the full Fock state \*S?), one finds 



B " (0) = y n \{N-n)\ (^-(0)^(0))* . (8) 

We consider first an approximation when Pij(t) can 
be found analytically. Specifically, we ignore the off- 
diagonal terms of the Hamiltonian (5) which correspond 
to the exchange of atoms between the states 0, 1 in the 
representation (4). Below it will be seen that this ap- 
proximation provides qualitatively correct results when 
these exchange terms are smaller than the term en in 
H n ,n (5). In future work we will show how this approx- 
imation can be extended for the case of larger values of 
the off-diagonal terms. Correspondingly, one has B n (t) — 
B n (0)exp(—iH ntn t). This results in Eq.(6) rewritten as 
poo = A r |C (0)| 2 '=const, p n = A|Ci(0)| 2 =const and 

Pw (t) =Npiexp(-i(e + ffl /2)t) _ 

■(l + ldWI^expH^-l))^ 1 , [J) 

with pi — C*(0)Co(0). Note that a closed form solu- 
tion similar to Eq.(9) can be obtained for a system with 
arbitrary numbers of levels in the same approximation 
neglecting the exchange effects. We will not discuss this 
case here. 

For N » 1 the short time behavior of Eq.(9) is gaus- 
sian, i.e. pio(i) = Apiexp((— ito\t — (£/r c ) 2 ), where we 
have introduced the effective Rabi frequency and the in- 
verse of the collapse time 

eoi=eii-£oo + |Ci(0)| 2 iV 5 i, - = ^AVn\ 9i \., (10) 

T~c V2 

respectively, and o(l/N) terms have been neglected. 
Note that eoi and r c are sensitive to the initial condi- 
tions given by pi. As a matter of fact, the dependence of 
eoi on the square of the initial amplitude (|Ci(0)| 2 « p 2 
for small |Ci(0)|) agrees with the observations of the 
parabolic dependence of the mode frequency on the re- 
sponse amplitude [&[§. Furthermore, r c in Eq.(10) is 
a nonanalytical function of p\ as well as the interac- 
tion constant g\ ~ Uq. This implies that the effect 
of the dephasing can not be derived by any sort of 
perturbation expansion starting from the noninteract- 
ing picture. Note also that the classical Eq.(7) in the 
same approximation gives p q' = poo, Pi* — Pn and 
p^\t > 0) = A r p 1 exp(— ie m t). Of course, this behavior 
is purely harmonic and demonstrates no signs of dephas- 
ing. At short times t « t c the MF solution Pi* (*) 



basically coincides with the quantum mechanical expres- 
sion (9). 

The solution (9) is sug gestive of another peculiar fea- 
ture - the revivals |^2| , pT| of the Rabi oscillations with the 
period t t = 2t: / ' g\ » t c . Below we will show that these 
main features of the analytical approximation which ne- 
glects the exchange terms remain qualitatively correct for 
the exact numerical solutions of the quantum as well as 
classical Eqs.(5), (7), respectively. 

We have diagonalized the Hamiltonian (5) numerically 
and then calculated the density matrix (6) for the initial 
condition (8). For N < 200 it is possible to do this di- 
rectly. For larger N the required calculation time grows 
significantly with N. Utilizing the fact that the initial 
condition (8) peaks around some value of n allowed us 
to reduce the size of the matrix (5). As a matter of 
fact, the distribution (8) is the square root of the bi- 
nomial distribution whose width is |Ci|yI — \Cq\ 2 >/N. 
This implies that for large N and small C\ the effective 
size of the Hamiltonian matrix grows as ~ \C\\\fN in- 
stead of as N. Accordingly, we reduced the size of the 
matrix (5) around the peak values of the distribution (8) 
and ran the program for several successively increasing 
sizes until the final solution did not change. The re- 
sult is shown in Fig.l for two different initial amplitudes. 
For comparison, we have also plotted the envelope of 
the classical solution of Eq.(7) found numerically as well. 
The Hartree ground state function i/jq and the function 
ipi are represented by the variational ansatz -0ooo( x ) = 
(M/7r) 3 / 4 wy 4 wy 2 exp(-Afw ± (a; 2 + y 2 )/2 - Muj z z 2 /2) 
|Ll| and the state ^iio(x) = (Mlu ± (x 2 + y 2 ) - l)^ooo(x), 
respectively. Here w z , w±_ stand for the variational fre- 
quencies |ll|| along and perpendicular to the axial axis 
z, respectively. These states were taken for given N. 
Accordingly, one finds for gi in Eqs.(5),(9), (10) g\ — 
—a\/Muj z uij_/^/2TT. The bare energy of the single parti- 
cle transition is £\\— £oo = uj±+ui 2 /uj±. The solutions for 
uj z . uj±_ were taken from Ref . |ll[] , with the frequency lu± 
at N = being lj = 2tt ■ 132Hz f|. When N < 300 such 
a choice corresponds to the lowest breathing mode, which 
for large N should transform into that found in Ref. ||. 
For N > 500, the use of the pair of states V'ooo, V'no f° r 
calculating the matrix elements in (4) becomes a poor 
approximation for the actual breathing mode j^] . There- 
fore in this paper we will limit our analysis to the case 
N = 300. 

We calculated the quantum mean square radius 
i? 2 (i) = (A/Too) -1 r ijPij(t) m tne xy-pl&ne of the 
condensate in units of Tq = 1/Muj±, where explic- 
i% ^oo = I d ^ooo( x2 + y 2 Hooo, r 2 x = J dxt/j* 110 (x 2 + 
y 2 )ipiw = 37-oQ and r 2 = / dxi/) 110 (a; 2 +?; 2 )V'ooo = ^oo; as 
a physical quantity responding to the Rabi oscillations. 
The MF approach gives R 2 (t) = l + 2(|Ci| 2 + Re(C 1 *C )). 
Obviously, in the Hartree ground state (Ci = 0) this 
quantity equals 1. Furthermore, exact calculations of R 2 
for the true many body ground state of the Hamilto- 
nian (5) for N < 500 give a value which is very close to 
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1 as well. For the initial condition (8), one finds that 
i? 2 (0) - 1 = 2pi + 2|Ci(0)| 2 for both the Hartree state 
and the actual state governed by the Schrodinger equa- 
tion (5). We will analyze the case of small disturbances of 
the ground state. This corresponds to \C\\ < 1/3 which 
allows one to take R 2 (0) -1» 2pi with a very good pre- 
cision. Therefore the quantity 2p\ in Eqs.(9),(10) plays 
the role of the initial amplitude of R 2 (t) for the small 
amplitudes. 




time (sec) 

FIG. 1. Collapses and revivals of oscillations of the mean 
square radius R 2 (t) of the condensate for N = 300: a) the 
initial amplitude 2 pi — 0.1; b) the initial amplitude 2pi = 0.6. 
Because of large time scale chosen, the rapid oscillations are 
seen as a black background. The dashed lines correspond to 
the envelope of the classical solutions 

Note that the amplitude of the classical solution (see 
the dashed line in Fig.l) stays constant and equal to 2pi. 
In contrast, the collapses and the revivals of the ampli- 
tude of the oscillations of R 2 (t) determined quantum me- 
chanically are clearly seen in qualitative accord with the 
analytical solution (9) . We have approximated the initial 
collapse stage of the numerical solution by a gaussian fit 
to the envelope and found the characteristic collapse time 
t c as a function of the initial amplitude 2p±. The result 
for 1/t c is represented in Fig. 2. As can be seen, this graph 
is very close to the straight line predicted by Eq.(10). Its 
slope is 7.0s _1 while the analytical formula (10) gives 
10s _1 . These values yield the shortest dephasing time 
(corresponding to p\ — 1/2 in Eqs.(9),(10)) as « 0.1s 
which is comparable to the relaxation times (« 150 ms 
H and sw 250 ms |J) reported experimentally. 

A direct confirmation of the proposed mechanism 
would be the experimental observation of the revivals. In 
this regard we note that the expected values of r r are of 
the order of seconds as given by the analytical estimates 
and the numerical solution (see Fig.l). This is compa- 
rable with the lifetime of the condensate in the traps 
. Furthermore, the normal component which intro- 
duces irreversible dissipation would inhibit the revivals 



unless its concentration is strongly suppressed. We will 
discuss this point in greater detail in future work. An 
estimate of the degree of the required suppression can 
be obtained from the following considerations. The rate 
1/t„ of the thermal dissipation should be proportional 
to the relative amount p of the normal component in 
the trap. In other words l/r« = p/tq, where To is the 
relaxation time above Tbec when all particles consti- 
tute the normal component (p = 1). This time is about 
t « 20ms Consequently, in order to observe the 

revival of the oscillations at times t > T r , the tempera- 
ture should be reduced so that the condition t„ >> r r 
or p « T Q /r r ~ 10~ 3 holds. 
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FIG. 2. The rate of the first collapse as a function of the 
initial amplitude 2pi for N = 300. The line is a guide for 
eyes. 

An indirect indication would be the observation of the 
relaxation of the oscillations described by the gaussian 
dependence, with the relaxation rate proportional to the 
initial response amplitude (for the small amplitudes) as 
predicted by Eq.(10) and Fig. 2. Furthermore, Eq.(10) 
predicts a specific dependence on A. Let us focus on 
this fact in detail. The coefficient gi in Eq.(10) is given 
by the effective volume V occupied by the condensate 
as gi ~ a/V ~ ap/N. For small A the density p ~ A 
yielding l/r c ~ VA in Eq.(10). For N > 200 — 500 the 
volume V expands as a function of A due to the inter 
particle interaction p| . Accordingly, the dependence of 
1/t c should deviate from ~ y/~N. If p were independent 
on A, one would have obtained l/r c ~ p/VN — > which 
implies no dephasing for large A. In the approximation 
we employed contributions of the higher order terms in 
p to g\ were neglected. Therefore, even though p grows 
with N (p ~ TV 2 / 5 |ll| for large A), the dephasing rate, 
as given by Eq.(10), is expected to vanish as a slow func- 
tion of A (1/t c ~ A -1 / 10 ) for large A. The higher or- 
der terms in p contributing to the dephasing rate could 
change this conclusion. In future work we will analyze 
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the [imit of large N by means of finding a better choice 



for the single particle state ipi for large N using the Bo- 
golubov approach. 

In conclusion, we have proposed a mechanism of inter- 
action induced dephasing of the collective excitations in 
the atomic traps. The calculated collapse time is compa- 
rable to the observed relaxation times for the collective 
modes. Two main features - gaussian relaxation with 
the time constant depending linearly on the initial am- 
plitude of the oscillations and depending on the total 
number of atoms - if observed experimentally, could ver- 
ify this mechanism. If the temperature could be reduced 
much below Tbec, the revivals of the collective oscilla- 
tions could be seen on times of several seconds. 
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